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In this paper we follow the analysis and protocols of recent experiments, combined with simple 
theory, to arrive at a physical understanding of quasi-condensation in two dimensional Fermi gases. 
We find that quasi-condensation mirrors Berezinskii-Kosterlitz-Thouless behavior in many ways, 
including the emergence of a strong zero momentum peak in the pair momentum distribution. 
Importantly, the disappearance of this quasi-condensate occurs at a reasonably well defined crossover 
temperature. The resulting phase diagram, pair momentum distribution, and algebraic power law 
decay are compatible with recent experiments throughout the continuum from BEC to BCS. 


Understanding two dimensional (2D) fermionic super¬ 
fluidity has a long history relating to the Mermin-Wagner 
theorem [T] and Berezinskii [5], Kosterlitz and Thouless 
(BKT) physics |3]- More recently it has been viewed 
as important for addressing the phase fluctuation pic¬ 
ture (and related pseudogap phenomena) associated with 
high-Tc superconductors [3]. Current interest in 2D 
bosonic superfluids in ultracold atomic gases has revealed 
a general consistency with the BKT transition [SH7] . For 
2D fermionic superconductors and superfluids, however, 
it should be emphasized that there is some historical 
controversy [5] (beginning with Kosterlitz and Thouless 
0) surrounding observable signatures and applicability 
of BKT physics. 

Thus recent reports linn] of a form of pair conden¬ 
sation in 2D fermionic gases are particularly exciting. 
These follow earlier work addressing the ground state m 
and the higher temperature regime, away from conden¬ 
sation [T^. These experiments pnni show that strong 
normal state pairing is an essential component of 2D 
Fermi superfluids, even in the BCS regime. In fact, much 
of the theory invoked to explain these experiments was 
based upon true Bose systems. A characteristic feature 
of 2D superfluidity at finite T is the presence of nar¬ 
row peaks in the momentum distribution of the pairs, 
without macroscopic occupation of the zero momentum 
state. Throughout the paper this will be our definition 
of “quasi-condensation.” This quasi-condensation in mo¬ 
mentum space is associated with algebraic decay of coher¬ 
ence in real space. Importantly, the BKT-related tran¬ 
sition temperature is manifested as a sudden change in 
slope of a normalized peak momentum distribution for 
pairs. 

In this paper we present a theory of a 2D Fermi gas 
near quantum degeneracy and show how it reproduces 
rather well the results of these recent experiments 0110] 
through an analysis of the phase diagram, the pair mo¬ 
mentum distribution and algebraic power laws. Given 
the ground breaking nature of the experiments, it is im¬ 
portant to have an accompanying theoretical study which 
follows exactly the same protocols without any adjust¬ 
ments or phenomenology. Our approach is to be dis¬ 


tinguished from other studies of 2D Fermi gases msi- 
[53|. In particular, those addressing BKT physics msi- 
uniiinj, use existing formulae [M] and determine the 
unknown parameters to obtain By contrast here 

we reverse the procedure and follow experimental pro¬ 
tocols to thereby provide a new formula, involving com¬ 
posite bosons, for the transition temperature associated 
with quasi-condensation. In the homogeneous case, this 
is analytically tractable and presented as Eq. ([^ below. 

Importantly, there is a rather abrupt crossover out of 
a quasi-condensed phase at a fairly well defined tempera¬ 
ture Tqc. In the BEC regime this matches earlier theoret¬ 
ical estimates of the BKT transition temperature which 
are based on different theoretical formalisms IlIIMSj. 
We find that Tqc varies continuously with scattering 
length and, in reasonable agreement with experiment 0, 
the transition appears at a slightly higher temperature 
for more BCS-like systems. We infer that the physics 
driving this quasi-condensation derives from implications 
of the Mermin-Wagner theorem; that is, from the inabil¬ 
ity to condense except at zero temperature. To minimize 
the free energy, the system remains quasi-condensed for 
a range of finite temperatures. Since we, as in Ref. [25], 
make no reference to vortices we cannot argue that our 
observations correspond strictly to a BKT scenario 0, 
but we can establish that our findings follow rather pre¬ 
cisely those of recent experiments. 

Background theory.— Theoretical studies of the 2D 
Eermi gas divide into two classes: those which build 
on or extend BCS mean-held theory jij [TBHTFl HSl US] j 
which is the largest class, and those (based on t-matrix 
schemes) which do not [T7IE(71 ITT] . Here we consider a 
t-matrix theory belonging to the hrst class. In the fol¬ 
lowing overview we omit technical details which can be 
found in the Supplemental Material and are extensively 
discussed elsewhere [551 US]. There, we also present a 
comparison with other theories. 

To describe the Eermi gas, we begin by introducing 
a pair propagator F(Q), representing a Green’s function 
for bosonic, or paired fermionic, degrees of freedom. Here 
we dehne the vector Q = (10, q), where iUl is a bosonic 
Matsubara frequency at temperature T and q is the pair 
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momentum. The pair propagator r(Q) is chosen so that 
r“^(0) = 0 at a temperature below a true 3D phase tran¬ 
sition temperature (where /ipair = 0 ) and, importantly, 
we impose the condition that this Thouless criterion re¬ 
produces the usual mean-field equation determining the 
pairing gap, A 0. We emphasize that A is a pairing 
gap and not an order parameter. In 2D, where Tc = 0, 
this equation at non-zero T is naturally generalized to 
r“^(0) oc /Ltpair- This serves to implement the reasonable 
assumption that we consider the normal phase for T > 0 
without phase coherence, but in the presence of a pairing 
gap so that A 7 ^ 0. 

A key component of the theory is the inclusion of fluc¬ 
tuations, or bosonic degrees of freedom. As we will show, 
fluctuations in 2D are necessarily unable to condense, 
thus guaranteeing that /ipair will never vanish for any 
T > 0. Because r((3) represents a pair propagator, it 
can be expanded at small Q into the generic form: 


r(Q) 


a 


-1 

0 


in riq -t- /ipair T i^Q 


( 1 ) 


and we associate Dq « /2Mb with a pair dispersion 
of mass Mb- Throughout we find that we can drop the 
small lifetime contribution yg. Note that the small Q 
form of the pair propagator is, up to a constant og, that 
of a Bose gas which has no direct inter-boson interac¬ 
tions, but in which the bosons interact indirectly via the 
fermionic medium. 

Performing the sum over bosonic Matsubara frequen¬ 
cies iQ gives the momentum distribution of bosons de¬ 
fined through ns (q) = aoI]inr((3) = 5 (flq - /tpair), 
where b{x) = {c^/^bT _ jg |;]qg Bose-Einstein distri¬ 
bution function. From here it is natural to define a boson 
number density ub through ub = X^q^s (‘l) = «oA^, 
where in the second equality we have associated the num¬ 
ber density with the pairing gap. This association is 
based on the self-energy and addressed in detail in the 
Supplemental Material, where we review the microscopic 
basis [551HH] of our theory. 

With this formalism we can now determine the un¬ 
knowns that appear in Eq. Q. We use the generalized 
gap equation r“^(0) = ag/tpair, along with the usual BCS 
equation for determining the fermionic chemical poten¬ 
tial, /i and the bosonic number equation to arrive at three 
coupled equations 
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= n, 


( 2 ) 
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(4) 


where n is the density of fermions. Here we de¬ 
fine the Fermi-Dirac distribution function / [x) = 


[e^lkeT _|_ single particle dispersion ^k = 

k'^ j2m — n for a fermion of mass m, momentum k and 
chemical potential fi, and the Bogoliubov dispersion with 
gap A which is given by Ak = have 

regularized the gap equation in Eq. ® by introducing 
a two particle bound state energy cb = h^/ma^D [23j . 
To match with experiment, we use a quasi-2D scattering 
length a 2 D parameterized through ln(fc 7 ?a 2 D)- We as¬ 
sume throughout that the transverse confinement is suffi¬ 
cient we can neglect corrections due to a finite transverse 
trapping length m- 

At T = 0, we can use the well known solution fi = 
ef — £ 5 / 2 , A = \J2tp(.B [23] along with /tpau = 0, where 
tp = TTh^n/m is the Fermi energy. At finite tempera¬ 
tures, Eq. ^ can be inverted exactly to give: 

/ipair = fcBTIn (5) 

where \b = \/2'Kh?'/MBkpT is the thermal wavelength 
for the bosonic pairs. The pair chemical potential there¬ 
fore crucially relies on the bosonic phase-space density 
"Db = ~ ^/T, which functions roughly as a proxy 

for inverse temperature. At low temperature, Vb ^ 1 
and we find /ipair/feB?" or that the chemi¬ 

cal potential is exponentially suppressed. On the other 
hand, at high temperatures <C 1 and /ipair ~ —TlnT 
which can be substantial. 

Analysis.— Since /ipair is finite and continuous at all 
non-zero temperatures there can be no true phase transi¬ 
tion m- Nevertheless there is a rather abrupt threshold 
from a moderately large to an exponentially small chemi¬ 
cal potential. We introduce a tolerance factor e which can 
ultimately be determined from the experimental proto¬ 
cols [5] , and which defines this threshold via the fugacity 
^ _ gAipair/fcsT^ boundary between the low and high 
temperature behaviors occurs when the slope of the fu¬ 
gacity with respect to phase-space density is of order e: 
dz{'DB)/d'DB ^ e. This will introduce a scale for an ef¬ 
fective (BKT-like) crossover temperature: Vb = ln(l/e). 

Importantly, there is a very weak (logarithmic) de¬ 
pendence on this tolerance factor which underlines the 
fact that the transition will be quite abrupt. More 
accessible experimentally |5] than /Xpair is the behav¬ 
ior of the zero pair momentum peak magnitude, called 
nB(q = 0) = nB(0). The magnitude of this peak is di¬ 
rectly related to the bosonic phase space density through 
tt-b (0) = — I. Thus we can rewrite the crossover 

constraint on the fugacity as dnB{0)/dVB ~ 1/e. In this 
way we find that a threshold in /tpair enters as a slightly 
rounded knee in 71^(0). Here we use this knee to deter¬ 
mine an effective quasi-condensation transition temper¬ 
ature, in much the same way as in experiment. In more 
direct comparisons with Ref. |9] we find that e is ap¬ 
proximately one percent. This corresponds to Vb ~ 4.6, 
which is close to the Monte Carlo result for the BKT 
transition of a true bosonic gas jSlj , which took a typical 
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FIG. 1. (a) Pair chemical potential, /ipair, as a function of temperature at a scattering length of ln(fcF 02 D) = ±0.5. At high 
temperatures, the pair chemical potential is large and negative, whereas it is exponentially suppressed as T —> 0. The inset 
shows the same quantity on a logarithmic scale to show the exponentially small but, non-zero, behavior, (b) Normalized 
momentum distribution at T/Tf = 0.07 and ln(fcFa2D) = —0.5. In the main plot, the red (solid) curve shows the trap-averaged 
pair momentum distribution in the high temperature limit. The black (dashed) curve is the trap-averaged interpolation from 
a corresponding Boltzmann distribution. The inset shows the same quantity from Ref. [^; the dashed line is a Boltzmann fit 
to the experimental data at large q. (c) Peaks in the normalized momentum density as functions of temperature at a range of 
scattering lengths across the BEC-BCS crossover. The black solid circles are Tqc determined from Eq. The insets on the 
right show an example of corresponding experimental data with ln(fcFa2D) = —0.5. 


value of Vb ~ 4.9 in Ref. nnj. 

Converting the transition condition to a form analo¬ 
gous to the widely applied [H [TiHTH] Kosterlitz-Nelson 
condition [24] . ksTc = ^Ps{Tc)k?/rn? for a phase stiff¬ 
ness Ps, we find: 


TT fi^nBiTqc) 
2.3 MB(rqc) ■ 


( 6 ) 


In the deep-BEC regime (where ub/Mb = njim) this 
yields kBTqc ~ which is similar to estimates in the 
literature given by = ^ep- Importantly, the 

present expression for Tq^ applies throughout the BCS- 
BEC crossover. Towards the BCS limit the number of 
bosons decreases, but this is compensated largely in the 
crossover temperature by the decrease in bosonic effective 
mass. These analytic arguments apply to a homogeneous 
system and, following experimental protocols, they relate 
to the characteristics of the non-condensed pairs. 

In order to better compare to experiment, we apply 
the local density approximation (LDA) to account for 
trap effects. We note the presence of a trap provides 
only a minor quantitative change to the general qualita¬ 
tive picture. To apply the LDA, we rewrite our equations 
using the transformations p ^ p (R) = po — 
and A —)■ A (R), where R is a local position, and is not 
to be confused with the conjugate to q. Note that po is 
still homogeneous in space (i.e., there is only one degree 
of freedom) but that A (R) is no longer homogeneous. 
We now have fk —t Ck (R) = k'^ /2m — p (R), and Ek —f 
Ek (R) = VCk (R) + (R)- Similarly, Mb ^ Mb (R), 

ub ^ ns (R), Go oq (R), etc., through these same 
substitutions. We also define a trap-integrated momen¬ 
tum distribution: ^^(k) = f ^^(k, R)d^R. For details 
of the LDA, including the parameters used, see the Sup¬ 
plemental Material. 


Comparison with experiment.— We now compare our 
theory with the recent experimental results in Refs. IS 
ITT)] , using our numerical results for the trapped case as 
“data” analogous to the experiment. In order to probe 
the momentum distribution of bosonic pairs at low tem¬ 
peratures, in Fig.[^a) we plot the pair chemical potential 
versus temperature for two values of ln(A:Fa 2 D) = ±0.5, 
along with an enlarged plot of \ppa,ir{T)\ which is pre¬ 
sented in the inset. The dots illustrate the crossover 
points associated with the threshold discussed earlier; 
here Pp^ir begins to appreciably deviate from zero, thus 
marking the transition out of the quasi-condensed state. 
Figure [^b) shows an example of the trap-integrated pair 
momentum distribution nB(q) at T < Tqc- The small 
chemical potential ^pair results in a sharply peaked dis¬ 
tribution hB(q) as q — 0. This behavior is similar to the 
results observed in experiment [S], as shown in the inset. 
Therefore this peak, or signature of quasi-condensation, 
emerges when the pair chemical potential becomes suffi¬ 
ciently small. 

To quantify an effective crossover transition tempera¬ 
ture, we focus on the ratio of the peak magnitude 11^(0) 
of this momentum distribution normalized to the peak 
number density in the center of the trap, no = n(R = 0), 
following the experimental protocol [5]. This is plotted 
for three different values of ln(/cFa 2 D) in Fig. [^c) with 
the dots indicating the knee (assuming a one percent tol¬ 
erance factor). This allows us to arrive at a BKT-like 
transition temperature in a trap as a function of scatter¬ 
ing length. The inset plots the experimental results for 
comparison. For our “data” the peak of the momentum 
distribution grows exponentially at low temperatures. In 
contrast, the experimental data does not grow quickly 
enough to distinguish between exponential growth and a 
sharp peak in this distribution. More generally, we note 
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(a) Phase Diagram 
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FIG. 2. (a) The phase diagram of a quasi-condensation 

crossover in a 2D fermi gas as a function of the scattering 
length ln(fc_Fa 2 D). The black curve in the phase diagram 
represents Tqc, and the colored shading represents the non- 
thermal fraction (see text and Fig. Bb))- The main structure 
can be understood from Eq. along with the pair mass (b) 
A/s(Tqc) and bosonic number density (c) UBiTqc) at Tqc. On 
the BEC side, Tqc limits to a constant value. A maximum of 
the pair mass near ln(fcFa 2 D) = 0 results in a drop in Tqc at 
the same point. In the BCS limit the transition temperature 
drops as the bosonic number density declines more quickly 
than the pair mass. All quantities are calculated from i? = 0 
data In the LDA. 


that the deviations between experiment and our theory 
consistently suggest that the absolute value of our chem¬ 
ical potential is too small. This, in turn, reflects the be¬ 
havior of the BCS-like gap equation which sets the scale 
for ^pair through Eq. ([2]). 

We next use this analysis to obtain the phase dia¬ 
gram as a function of interaction strengths investigated 
in experiment. The results are shown in Fig. Ba) plot¬ 
ted against ln(A:i;’a 2 D)- We can associate a BCS-like 
phase with positive fermionic /r which appears when 
ln(fc^a 2 D) > 0. Figures^b) andj^c) indicate the numer¬ 
ator {ub) and denominator {Mb) components of Tqc as 
shown in Eq. § as a function of scattering length. The 
color coding indicates the non-thermal fraction which is 
found from Fig. as the area between the momen¬ 
tum distribution (solid curve) and its high temperature 
asymptote (dashed line). On the BEC side, we find that 
Tqc saturates. As the scattering length is increased, the 
transition temperature begins to drop before rising to 
a local maximum and then falling off in the deep BCS 
regime. While the values are rather similar, this non¬ 
monotonic behavior is not as directly seen in experiment, 
although it is suggested in their plots of the non-thermal 
fraction. It should also be noted that our theory for the 
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(b) Experimental First-Order Correlation Function 



FIG. 3. Correlation function, (/i(r), at large distances found 
through a Fourier transform of the trap-averaged momentum- 
density in the LDA approximation. Both the theoretical and 
experimental results are presented on a log-log scale. In the 
top panel, we consider three different temperatures at a fixed 
scattering length ln(fe_Fa 2 D) = —0.5, where Tqc = 0.094 Tf. 
From top to bottom, T/Tp = 0.06, T/Tp = 0.09, and 
T/Tp — 0.12. For T < Tqc, we ht to a power law expres¬ 
sion (dashed lines). The dashed curve for T/Tp = 0.12 is an 
exponential fit. The bottom panel is the experimental com¬ 
parison. The temperatures for the theory curves are close to 
those considered in experiment. 


transition temperature is valid across the entire BEC- 
BCS spectrum, rather than as two endpoint cases as of¬ 
ten studied [HI [ 20 ]. 

Finally, in Fig. we present the correlation function 
(/i(r) determined from the Fourier transform of the trap 
integrated momentum distribution hB(q), again follow¬ 
ing the experimental protocol unj. We fit to a power 
law for an intermediate range of r corresponding roughly 
to that used in the experimental data [10]. It should 
be noted, however, that the power law regime appears 
slightly more extended in experiment than in theory. 
With our analytic insight we believe there may be bet¬ 
ter fits to our “data” by considering the structure of 
the momentum distribution at q —)■ 0 [33]. Neverthe¬ 
less following experiment, we find a reasonable fit to a 
power law in this range at low temperature, gi(r) ^ l/r’^, 
and a crossover to an exponential fit at higher tempera¬ 
tures, gi(r) ^ Our power laws lie in the range of 

0.8 < r] < 1.45. These values are close to the power laws 
observed in the experiment of 0.6 < 77 < 1.4, which ap¬ 
pear universal near the expected BKT transition temper¬ 
ature, yet are far from the predicted scaling of ry < 1/4. 

Conclusions.— The favorable comparisons between 
theory and experiment in Figs. B 1 provide helpful in¬ 
sights into the behavior of 2D Fermi gases. Central to 
our picture, is the relation between the zero momen¬ 
tum peak in the pair distribution function and the small 
pair chemical potential yxpair- As consistent with the 
Mermin-Wagner theorem, /ipair is shown to never van- 
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ish except at zero temperature. We argue that it is this 
inability to fully condense which ultimately drives quasi¬ 
condensation. Importantly, with increasing tempera¬ 
ture there is a rather abrupt transition from this quasi- 
condensed phase which, following experimental protocols 
[siiin] is reflective of BKT physics. 

Our approach should be contrasted with contributions 
to the theoretical literature which address BKT physics 
muSHTllEo], by solving for the phase stiffness, etc. pa¬ 
rameters that appear in the usual formulae [Ml US] for the 
BKT transition temperature. Here we reverse the logic in 
order to present an alternative viewpoint. We follow the 
experimental procedure to thereby provide a new formula 
(see Eq. in the absence of a trap) for the transition 
temperature associated with quasi-condensation. This 
can be expressed in terms of the bosonic phase space 
density, Vb, which we find, as in other literature |32j . 
is around four; as expected m, "Db exhibits no sharp 
transition. We stress that our expression is associated 
with composite bosons whose mass and number density 
vary significantly and continuously from BCS to BEC. 
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I. OVERVIEW OF THE PRESENT PSEUDOGAP FORMALISM 


The theoretical formalism in this paper is a t-matrix theory chosen to reproduce the physics of a generalized BCS 
phase. This approach establishes the nature of non-condensed pairs of fermions, which will self-consistently generalize 
BCS theory. We confine ourselves to two dimensions where, as we have seen in the main text, there is no condensate. 
All fermion pairs are thus non-condensed. The BCS dispersion relation can be written as = ■\/Ck + j where we 
write the gap parameter as A and make it clear that this represents a pairing gap and not an order parameter. In all 
t-matrix theories the t-matrix, r((5), is written as 


r(Q) 


9 

1 -I- gX{Q) 


(SI) 


Each t-matrix theory differs only in the choice of the pair susceptibility, X{Q). 

To make the appropriate choice, here we note that the mean-field gap equation can be expressed as T”^ (0) = 0, 
where 




f{E^) - /(gk) 


(S2) 


Note that this is associated with a pair of Green’s functions, one of which is dressed and one of which is bare. We 
now write a natural generalization of r(0) to r((5). This leads to 


r-'(0) = E 


r 1 - /(Au) - /(a-q) 2 _ /(Ak) - /(a-q) 

- Ak -f Ck-q — — z0+ ^ E\g_ — ^k-q + + *0+ 


(S3) 


which corresponds to xiQ) = J2k E'{K)Go{Q — K) where Go(Ar) is a non-interacting Green’s function, and the 
full Green’s function G{K) includes a self-energy T,{K) through G~^{K) = Gq^{K) — E,{K). The specific form of 
Eq. (S3) requires a BGS-like self-energy of the form T,{K) = —A^Go(—A). Note that the general GGq form of x can 
be derived using the equations of motion for the four-point correlation function in the ladder diagram approximation 
and is known to be the only t-matrix consistent with the BGS gap equation [SI, S2]. 

We now show that the above t-matrix theory is consistent with a BGS-like pairing gap in the limit that p.pair is 
sufficiently small. As in all t-matrix theories, the fermions acquire a self-energy which is given by 

S(A) = ^r(Q)Go(Q-A). (S4) 

Q 


If the pairs are close to condensation (with T ^(0) « 0, so that r(0) « oo) then we can approximate the above 
equation by 


E(A) 



Go(-A) 


-A^Go{-K), 


(S5) 


where in the second equality we have defined A^ = — r((5) to show that this approximation is consistent with a 

BGS-like self-energy. Note that A^ enters only as a pairing gap in the Bogoliubov spectrum Ak, and does not need 
to represent an order parameter in any way. Thus, we have now established a self-consistent t-matrix formalism that 
reproduces the BGS gap equation in the limit that = 0. 

At small Q, the t-matrix can be expanded in the form 


r(Q) 


a 


-1 

0 


^ ilq + /Xpair + *7 q’ 


(S6) 
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where the pair chemical potential is ^pair = ^ (s ^ + x(0)) and we introduce the bosonic effective mass from a 

small momentum pair dispersion flq « /2Mb- The specific expressions for Mb and 7 q are complicated and are the 

same as those in Refs. [SI, S2] up to a replacement of a 3D sum with a 2D sum. The constant gq = {n — no)/2A^, 
where n = 2 G{K) and Uq = 2 'YIk Go{K). From our definition of A^, after summation over bosonic Matsubara 
frequencies 

A2 = -^r(Q) = -^&(Dq-^pair), (S7) 

Q q 


where b(x) = i jg Bose-Einstein distribution function as dehned in the main text. This allows us to 

identify 


UB = aoA^ = -(n-no), 


(S8) 


where ub represents the total number of non-condensed bosons. 

In summary, the essential approximation in this scheme is the assumption that /Xpair is small, so that Eq. (S5) is 
valid. As a consequence, the fermionic self-energy is of the BCS-type, and the number equation is essentially that of 
BCS theory with the pairing gap A. 

These arguments lead to the three coupled equations in Eqs. (2—4) of the main text. However, we can use the 
definition of /Xpair in the text to reduce the system of equations to be solved for to: 


E 


1 - 2/ (Ek) 

2Ak 


1 

2ek + es 


E 1 


ill 

Ak 


(l-2/(Ak)) 


= aorin , 


= n. 


(S9) 

(SIO) 


This is now two equations in two unknowns. Here is the binding energy defined in the main text in terms of the 
scattering length and Xb = ^/2TrWjMBkBT is the thermal wavelength for the bosonic pairs. 


II. TRAP EFFECTS 


In the presence of a trap, the homogeneous mean-field equations are no longer valid. However, some trap effects can 
be included if we rewrite our equations using the transformations /i —>■ /i (R) = Ho — and A —> A (R), where 

R is a local position, and is not to be confused with the conjugate to q. Note that ho is still homogeneous in space 
(i.e., there is only one degree of freedom) but that A (R) is no longer homogeneous. We now have ^k ?k (R) = 
k"^/2m — H (R), and Ek —>■ (R) = ■\/^^~(R)~+~A2~(R). Similarly, Mb —>■ Mb (R), ub — > u-b (R), oq oq (R), etc. 

through these same substitutions. 

Using these substitutions, known as the local density approximation (EDA), the mean-field equations now become 


E 


1-2/ (Ak(A)) 
2Ak (R) 


1 

2ek -l- eb 


Ub (R) 
^(R) 


Tin 


^ 


27r 


dR R 


El 


^^^{l-2fiE^{R))) 


= N, 


(Sll) 


(S12) 


where N is the total number of particles, and we have used cylindrical symmetry. These equations are discretized 
on a grid of M points in R space. Equation (Sll) is split into M equations, whereas Eq. (S12) is a single equation. 
Therefore, we have M + 1 total equations. Furthermore, since we have discretized A [R] into M points with fio ^ 
single parameter, we have M -\- 1 degrees of freedom. Thus, there are an equal number of constraints as there are 
degrees of freedom. 

Within the EDA, the solutions to the mean-field gap equations at T=0 can be written analytically: 


A [R) = \/2eB^F 






1 - 



es 

T' 


(S13) 


H {R) = ep 


2 


(S14) 
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Here the Thomas-Fermi radius Rtf is defined by muj'^ p /2 = ep- 

In the LDA, the Fermi energy is defined from the T = 0 and A = 0 limit which, including the factor of two for 
spin degeneracy, gives ep = '/Nhw. From the LDA chemical potential /i(0) = ep, and we can connect the central 
density, n(0), to the Fermi energy in the usual way: ep = 1i^kp/2m, where = 27rn(0) dehnes the Fermi momentum 
hkp. Throughout our numerical calculations we use N = 200^ particles, consistent with the experimental value of 
N = 40,000 total atoms. This relation then implies tiujjeF = 1/200. 

In both the theory and the experiment, the scattering parameter is input as ln(fci^’a 2 D)- Note that the experiment 
estimates the Fermi energy through the central density at finite temperature leading to a slightly different value of 
ln(fci^’a 2 D) for an equivalent 020 - To show this deviation is not significant, we define r = ^/Nkw/{TTh‘^n{0)/'m), i.e., 
the ratio of the Fermi energy as defined from the trapping frequency and atom number, to the Fermi energy as defined 
from the central density. Using the numbers from Ref. [S3], we find r is in a range of r ~ 0.25 to r ~ 0.9 in the 
BEG and BCS limit respectively. Thus there may be a weak difference in comparing our scattering length parameter, 
ln(fci^’a 2 D), relative to the experiment [S3, S4]. We do not expect this correction to significantly modify our analysis 
or comparison with experiment. 


III. COMMENTS ON OTHER THEORETICAL APPROACHES 

Alternative theories of 2D superfluids at finite temperature (which treat BKT physics to varying degrees) fall 
into two main classes: (1) Theories which have a BCS basis and (using path integral techniques [S5-S9]) focus 
on determining the phase stiffness parameter, ps- This is then substituted into the BKT transition temperature 
expression. (2) Those which adopt one of two t-matrix approaches [S10-S13] neither of which, in contrast to the 
present paper, is based on BCS theory. 

The first class of theories [S6-S9], as noted in Ref. [S7], neglects non-condensed pair degrees of freedom [S14], 
in contrast to the second class [S10-S13], but in the first class, BKT physics plays a more prominent role. All 
homogeneous theories in 2D must satisfy the Mermin-Wagner theorem. For the first class of theories this theorem 
can be directly imposed by performing an average over phase fluctuations. As for the second class, we quote from 
Ref. [S7], that the Mermin-Wagner theorem “zs an excellent check on whether one has a good solution to the t-matrix 
equations: the correet solution should give Tc = 0 in 2D.” We have verified this consistency in the present t-matrix 
formalism. 

In addition, there are theories which introduce inhomogeneity, trap and/or finite size effects [S5, SI5] which thereby 
enable the calculation of a finite transition temperature often also associated with the BKT phase and there is also 
an extensive literature on the ground state [SI6, SI7]. 

As for alternative t-matrix theories, in Ref. [S13] a Baym-Luttinger-Ward approach is applied to study a 2D 
homogeneous Fermi gas. In the Baym-Luttinger-Ward formalism the self-energy arises as a functional derivative of 
the Luttinger-Ward functional with respect to the full Green’s function. The advantage of this formalism is that 
it is manifestly conserving and self-consistent. However, if one truncates this infinite system of diagrams, then the 
theory is no longer conserving. Thus, in 3D for example, satisfying the Thouless criterion requires a redefinition of 
the coupling constant parameters [SI8]. In 2D, where there is no Thouless criterion, one must regularize the theory 
to ensure the Mermin-Wagner theorem is satisfied [SI9]. Reference [S13] applies a condition that F“^(0) oc 1/N 7 ^ 0, 
which appears in the BEG literature [S20], and sets N = 500. Due to numerical difficulties they are only able to 
address this BKT transition on the BGS side of resonance. 

Recent work [Sll] based on a GqGq t-matrix formalism has presented results for the pair momentum distribution 
which appear to have some overlap with the present paper and the experimental data in Ref. [S3]. However, in 
contrast to both this paper and experiment. Ref. [Sll] does not include a trap, which makes a comparison between 
theory and experiment difficult. Indeed, in Ref. [S12] the authors examined the effects of a trap within the same 
theory and found a very high superfluid transition temperature (not related to BKT physics). 
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